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ABSTRACT 


We present high-resolution (R ~ 18,000), high signal-to-noise 2 pm spectra of two luminous 
X-ray flaring Class I protostars in the p Ophiuchi cloud acquired with the NIRSPEC spectrograp 
of the Keck II telescope. We present the first spectrum of a highly veiled, strongly accreting 
protostar which shows photospheric absorption features and demonstrates the stellar nature of 
its central core. We find the spectrum of the luminous ( L bol = 10 L©) protostellar source, 
YLW 15 to be stellar-like with numerous atomic and molecular absorption features, indicative 
of a K5 IV/ V spectral type and a continuum veiling r k = 3.0. Its derived stellar luminosity 
(3 L©) and stellar radius (3.1 Rq) are consistent with those of a 0.5 M© pre-main-sequence star. 
However, 70% of its bolometric luminosity is due to mass accretion, whose rate we estimate to be 
1 7 x 10' 6 M© yr- 1 We determine that excess infrared emission produced by the circumstellar 
accretion disk, the inner infalling envelope, and accretion shocks at the surface of the stellar core 
of YLW 15 all contribute significantly to its near-IR continuum veiling. Its rotational velocity v 
sin i = 50 km s~ l is comparable to those of flat-spectrum protostars but considerably higher than 
those of classical T Tauri stars in the p Oph cloud. The protostar may be magnetically coupled 
to its circumstellar disk at a radius of 2 - 3 R,. It is also plausible that this protostar can shed 
over half its angular momentum and evolve into a more slowly rotating classical T Tauri star 
by remaining coupled to its circumstellar disk (at increasing radius) as its accretion rate drops 
by an order of magnitude during the rapid transition between the Class I and Class II phases of 

evolution. . . 

The spectrum of WL 6 does not show any photospheric absorption features, and we estimate 

that its continuum veiling is r k > 4.6. Its low luminosity (2 L©) and high veiling dictate that 
its central protostar is very low mass, M ~ 0.1 M 0 . We also evaluate multi-epoch X-ray data 
along with these spectra and conclude that the X-ray variabilities of these sources are not directly 
related to their protostellar rotation velocities. 
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1. Introduction 

The discovery of accreting protostars has been 
a triumph for understanding low-mass star forma- 
tion. Infrared and sub-millimeter data have re- 
vealed that these objects are deeply embedded , 
cold, and luminous. The spectral energy distri- 
butions (SEDs) of such Class I protostars grad- 
ually rise from the near-IR, peak in the far-IR 
(A ~ 100/rm), and fall off in the sub-mm and mm 
wavelength regimes (e.g., Adams, Lada, & Shu 
1987). Observations and models have indicated 
that protostars consist of three primary compo- 
nents: a massive infalling envelope, a circumstellar 
accretion disk, and an embryonic stellar core. 

With a size of ~ 10 4 AU, the infalling envelope 
is the largest component of a protostellar object. 
It extends to within a few AU of the embryonic 
stellar core and absorbs and reprocesses all the vis- 
ible (and UV) as well as much of the near-infrared 
radiation emitted by the central protostar. The 
infalling envelope contributes the bulk of the ob- 
served far-infrared to millimeter-wave luminosity 
and a significant portion of the detected near- and 
mid-infrared emission. 

The circumstellar accretion disk can be as much 
as a few hundred AU in size and it contributes a 
significant portion of the observed mid- and near- 
infrared luminosity of the protostar. During the 
protostellar phase of evolution, the embryonic stel- 
lar core at its center grows primarily through ac- 
cretion of material from the circumstellar disk, 
which itself acquires its mass directly from the in- 
falling envelope. 

The embryonic stellar core at the center of the 
protostar is thought to be a stellar-like object a 
few solar radii in extent (e.g., Stahler, Shu, & 
Taam 1980). It is the source of the bulk of the 
protostellar luminosity which is generated by both 
a strong accretion shock at the stellar surface and 
the radiation of internal stellar energy. This inter- 
nal energy is provided by its self-gravity as well as 
the release of gravitational potential energy during 
the protostar’s assembly from accreting and in- 
falling material. The core’s accretion-dominated 
luminosity is radiated primarily at UV and vi- 
sual wavelengths, but most, if not all of this ra- 
diation is absorbed by dust in the infalling enve- 
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lope which renders the central protostar invisible 
at these wavelengths. As a result, little is known 
about the nature of the embryonic stellar cores at 
the hearts of protostellar objects. 

The observational technique of high-resolution 
near-IR. spectroscopy has the potential to reveal 
the nature of these heavily embedded central ob- 
jects. Indeed, we have recently been able to 
demonstrate the stellar nature of the central ob- 
jects in flat-spectrum protostars and thus discern 
their spectral types, effective temperatures, rota- 
tion velocities, luminosity classes, and masses us- 
ing this technique (Greene & Lada 1997, hereafter 
Paper I). These objects are in a later evolutionary 
state than Class I protostars; they are still sur- 
rounded by substantial circumstellar disks and en- 
velopes but their accretion luminosities (and rates) 
have decreased below Class I levels. However, our 
previous high-resolution IR spectroscopic study of 
Class I young stellar objects (YSOs) did not re- 
veal any spectral absorption features in these pro- 
tostars (Greene k Lada 2000, hereafter Paper II). 
Luhman k Rieke (1999) did detect absorption fea- 
tures in the moderate-resolution near-IR spectrum 
of at least one Class I YSO (YLW 16A / IRS 44), 
but the moderate continuum veiling of that object 
( r k = 1) suggests that it is intrinsically similar to a 
flat-spectrum YSO. Also, Kenyon et al. (1998) de- 
tected TiO absorptions in moderate-resolution op- 
tical spectra of three Class I YSOs in the Tau-Aur 
dark clouds. However, these are all low luminosity 
objects, L, ~ Lt 0 i < 0.5 Lq, indicating relatively 
low levels of accretion. Thus there is very little or 
no existing information concerning the nature and 
basic physical characteristics of the central stel- 
lar cores of protostars undergoing significant mass 
accretion. 

Observational determination of the effective 
temperatures and radii (or surface gravities) of 
these deeply embedded central objects would pro- 
vide critical constraints for theoretical protostellar 
models which predict the values of these parame- 
ters as a protostar evolves (Stahler, Shu, k Taam 
1980). Empirical knowledge concerning the ini- 
tial angular momentum of these objects as well as 
the evolution of angular momentum in the earliest 
phases of young stellar evolution is also highly de- 
sirable to test protostellar theory. Many classical 
T Tauri stars (CTTSs) embedded in the p Oph 
cloud and elsewhere are known to rotate slowly, 
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tr sin i < 20 km s -1 . However, we have found that 
flat-spectrum protostars in this cloud have similar 
radii but rotate much more quickly, v sin i ~ 40 
km s -1 (Paper I). Do heavily-accreting Class I 
protostars in Ophiuchus rotate as fast or faster? 
If so, how is stellar angular momentum reduced so 
quickly by the classical T Tauri pre-main-sequence 
(PMS) evolutionary phase? 

We have undertaken a new, very sensitive near- 
IR spectroscopic study of Class I protostars with 
the 10 m Keck II telescope in order to address 
these issues. The Keck observations are consid- 
erably more sensitive than our previous deep in- 
frared spectral survey made with the 3-m NASA 
IRTF. In this paper we report our findings for two 
of these objects in the p Oph dark cloud, YLW 15 
(IRS 43) and WL 6. These are 2 of the 3 p Oph 
Class I protostars detected in hard X-rays by the 
ASCA satellite, and both showed strong, variable, 
and flaring emission (Kamata et al. 1997). We 
describe our data acquisition and reduction in §2, 
analyze the spectra in §3, discuss the results in §4, 
and summarize our conclusions in §5. 

2. Observations and Data Reduction 

Near-IR spectra of the p Oph protostar YLW 
15 were acquired on 2001 July 7 UT, and the pro- 
tostar WL 6 was observed on 2000 May 30 UT. 
MK spectral standards were also observed 2000 
May 29 - 30 UT. All data were acquired with the 
10-m Keck II telescope on Mauna Kea, Hawaii, 
using the NIRSPEC facility multi-order cryogenic 
echelle spectrograph (McLean et al. 1998). Spec- 
tra were acquired with a 0758 (4 pixel) wide slit, 
providing spectroscopic resolution R = X/SX — 
18,000 (16.7 km s -1 ). The plate scale was 0720 
pixel -1 along the 12" slit length, and the see- 
ing was 075 - 076. The NIRSPEC gratings were 
oriented to allow both the 2.21 Na and the 
2.30 ^m CO band head regions to fall onto the 
instrument’s 1024 x 1024 pixel InSb detector ar- 
ray, and its NIRSPEC-7 blocking filter was used. 
NIRSPEC was configured to acquire simultane- 
ously multiple cross-dispersed echelle orders 31 
- 36 (2.08 - 2.45 pm, non-continuous) in 2000 
May and orders 32 - 38 (1.97 - 2.38 pm, non- 
continuous) on 2001 July 7. Each order had an 
observed spectral range A A ~ A/67 (At; ~ 4450 
km s -1 ). 


The internal instrument slit rotator was used to 
maintain a fixed position angle on the sky when 
observing WL 6 and the MK standards on 2000 
May 29 - 30. The slit was held physically sta- 
tionary and thus allowed to rotate on the sky (as 
the non-equatorially-mounted telescope tracked) 
when observing on 2001 July 7 (YLW 15). Data 
were acquired in pairs of exposures of durations 
from less than 1 second (for giant MK standards) 
to up to 600 s (for the YSOs) each, with the tele- 
scope nodded 6" along the slit between frames 
so that object spectra were acquired in all expo- 
sures. WL 6 was observed for a total of 30.0 min- 
utes, while YLW 15 was observed for 73.3 minutes. 
Early-type (B7 - A2) dwarfs were observed for tel- 
luric correction of the MK standard stellar spec- 
tra. HR 6054 (B7 IV) and HR 6070 (A0 V) were 
observed for telluric correction of the WL 6 and 
YLW 15 spectra, respectively. The telescope was 
automatically guided with frequent images from 
the NIRSPEC internal SCAM IR camera during 
all exposures of more than several seconds dura- 
tion. Spectra of the internal NIRSPEC continuum 
lamp were taken for flat fields, and exposures of 
the its Ar, Ne, Kr, and Xe lamps were used for 
wavelength calibrations. 

All data were reduced with IRAF. First, object 
and sky frames were differenced and then divided 
by flat fields. Next, bad pixels were fixed via in- 
terpolation, and spectra were extracted with the 
APALL task. Extracted spectra were typically 4 
pixels (076) wide along the slit (spatial) direction 
at their half-intensity points. Spectra were wave- 
length calibrated using low-order fits to lines in the 
arc lamp exposures, and spectra at each slit po- 
sition of each object were co-added. Instrumental 
and atmospheric features were removed by divid- 
ing wavelength-calibrated object spectra by spec- 
tra of early-type stars observed at similar airmass 
at each slit position. Combined spectra were pro- 
duced by summing the spectra of both slit posi- 
tions for each object and then shifting the spectra 
in wavelength to correct for radial velocities. The 
standard spectra were shifted so that their lines 
matched laboratory wavelengths, while the YSO 
spectra were corrected for all Earth and solar sys- 
tem radial motions relative to the local standard of 
rest. The YSO spectra were not further corrected 
for the radial velocity of the p Oph cloud. Fi- 
nally, YSO spectra were multiplied by an artificial 


3 


spectrum of a T~10,000K black body in order to 
restore their true continuum slopes (although this 
was a small correction of under 10% per order). 

3. Data Analysis and Results 
3.1. Spectra 

The spectra of several MK standard stars are 
shown in Figure 1. The three NIRSPEC orders 
with the most numerous absorption features are 
displayed. These features are similar to ones seen 
by Wallace & Hinkle (1996) in stars of similar 
spectral type, although we present a finer grid 
of spectra in the K2V - M2V range. Spectra of 
these stars are dominated by strong features of 
Si (2.0923 pm), Mg (2.1066 pm), and A1 (2.1099 
pm) in the shortest wavelength order, with ratios 
Al/Mg and Al/Si being very sensitive to spectral 
type. Some of the strongest lines in the next or- 
der are Na (2.2062, 2.2090 pm), Sc (2.2058, 2.2071 
pm), and Ti (primarily 2.2217 and 2.2240 pm). 
The Ti lines are very good indicators of spectral 
type in this range, growing stronger from earlier 
to later types. The 2.2069 pm Si line is also a 
good indicator of spectral type, weakening from 
early to later types. Therefore the ratio of Ti/Si 
is a particularly good spectral type diagnostic in 
this order. The strongest features in the longest 
wavelength order are Mg (2.2814 pm) and the CO 
v — 0 — 2 band head and rotation-vibration lines 
(2.2935 pm and long-wards). The Mg line becomes 
significantly weaker in later spectral types. All 
transitions identified above and hereafter (and in 
Fig. 1) occur in neutral species. 

The two giant spectra (K3 III and M2 III) 
shown in Figure 1 have many of the same lines 
as the dwarfs, but there are significant differences. 
The Si and Mg lines in the shortest order do 
not appear to change monotonically with spectral 
type, but the lines in the next order are much 
more telling. The Sc (2.2058, 2.2071 pm) lines 
are very strong (absolutely and relative to Na) in 
the M2III giant, while the Ti lines are very strong 
in both giants. The Na lines (2.2062 and 2.2090 
pm) are less broad in the giants because they are 
formed under conditions of lower pressure than 
in dwarfs. The Mg (2.2814 pm) line is also less 
broad in the giants and does not vary monotoni- 
cally with spectral type. The 2.2935 pm CO band 
head and nearby rotation-vibration lines are very 


strong in the giants. Numerous other low-strength 
absorptions are seen in all three orders of the gi- 
ant star spectra, indicating the presence of various 
molecules (see also Wallace & Hinkle 1996). 

Three NIRSPEC spectral orders of the YSOs 
WL 6 and YLW 15 are shown in Figure 2. The WL 
6 data have signal- to-noise of 100 - 170 at most 
wavelengths, and the YLW 15 data have signal- 
to-noise 240 - 300. Both objects have strong con- 
tinua rising towards longer wavelengths. YLW 15 
has numerous absorption features which are weak 
and moderately broad, while any absorption fea- 
tures in the spectra of WL 6 are either weaker or 
non-existent. Both sources exhibit strong v = 1 
- 0 S(0) H 2 emission at 2.2233 pm. Interestingly, 
the strength of this emission appears to be corre- 
lated with infrared excess and veiling in these ob- 
jects. WL 6 has weaker (absent) absorption lines, 
a stronger red continuum, and stronger H 2 emis- 
sion than YLW 15. 

3.2. Veiling and Rotation Analysis 

The spectra of YLW 15 and WL 6 are com- 
pared to those of the K4 dwarf HR 5568 in Fig- 
ure 3. In this figure, the spectra of the YSOs are 
processed beyond that shown in Figure 2 in order 
to facilitate this comparison. First, their contin- 
uum slopes were removed by normalization. Next, 
the relative strengths of their absorption lines were 
magnified by subtracting off a constant continuum 
value of 0.75, 75% of the normalized continuum. 
This specific constant was chosen so that the ab- 
sorption lines would have nearly the same equiv- 
alent widths in both YLW 15 and HR 5568. The 
spectrum of HR 5568 was modified by convolv- 
ing it with a limb-darkened stellar rotation profile 
with velocity v sin i = 50 km s~ x . This was re- 
quired to match the velocity widths of its absorp- 
tion features to those of YLW 15. 

Figure 3 shows that the spectrum of YLW 15 
is very similar to that of HR 5568. Nearly all 
absorption features have nearly identical relative 
strengths. Comparing relative line strengths is 
most useful because ratios (but not absolute equiv- 
alent widths) are insensitive to continuum veiling. 
For example, the ratios of their Mg (2.1066 pm) 
and A1 (2.1099 pm) equivalent widths are very 
similar, as are the ratios of their lines in the Naand 
Sc region (2.2053 - 2.2096 pm). The lines in the 
Na and Sc region also have the same strengths rel- 
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ative to the Fe and Sc lines in the 2.2263 - 2.2273 
pm region. However, the spectra of the two ob- 
jects do not match exactly. The 2.2280 and 2.2317 
pm lines of Ti are stronger relative to the 2.2053 

- 2.2096 pm Na, Sc, and Si lines in YLW 15 than 
in HR 5568. Likewise, the 2.2814 pm Mg line is 
weaker relative to the CO band head (or the 2.2053 

- 2.2096 pm Na, Sc, and Si lines) in YLW 15 than 
in HR 5568. Both of these differences suggest that 
YLW 15 is of slightly later spectral type than HR 
5568. It is clear from Figure 1 that YLW 15 is not 
as late as G1 338B (K7 V) and that its spectrum 
is much more dwarf-like than giant-like (see also 
the discussion of giant spectra in §3.1). We assign 
the spectral type K5 IV/V to YLW 15. Differ- 
entiating between luminosity class IV and V will 
probably require comparisons to synthetic model 
spectra given the rarity of well-classified luminos- 
ity class IV stars. 

The continuum veiling of YLW 15 can be mea- 
sured now that its spectral type has been de- 
termined. Continuum veiling at wavelength A 
is defined as the value of continuum excess ex- 
pressed as a fraction of the stellar photospheric 
flux, r\ — F\ :ex /F\.. We have estimated the Ri- 
band continuum veiling of YLW 15, r kl by com- 
paring the equivalent width of its K-band spectral 
features to those of MK standards. The above con- 
tinuum veiling equation can be applied to K-band 
wavelengths and written in terms of the intrinsic 
(EW) and observed {EW) equivalent widths of 
spectral features there, EW - EW/ (l + r *). 

We have evaluated r k of YLW 15 by compar- 
ing the measured equivalent widths of its 2.2053 

- 2.2096 pm Na, Sc, and Si lines (from its spec- 
trum in Figure 2) to those of the K4 V star HR 
5568 and the K7 V star G1 338B. The measured 
equivalent widths of these features are 0.8, 2.8, and 
3.5 A for similarly-reduced spectra of YLW 15, 
HR 5568, and G1 338B, respectively. The intrinsic 
equivalent widths of these features is expected to 
be approximately 3.2 A for YLW 15, intermedi- 
ate between HR 5568 and G1 338B. Substituting 
these numbers into the above equation yields r k 
= 3.0. We estimate that this value is accurate to 
about 25% given the uncertainties in the measured 
equivalent widths and the spectral type determi- 
nation of YLW 15. This veiling implies that 75% 
of the K-band flux of YLW 15 arises from (cir- 
cumstellar) continuum excess, and 25% of its flux 


arises from its K5 stellar photosphere. Therefore 
the stellar spectrum of this object can be approx- 
imated by subtracting 75% of the continuum level 
at K-band wavelengths. This - and normaliza- 
tion by its continuum - was done to produce the 
spectrum of YLW 15 shown in Fig. 3. 

The derived veiling r k = 3.0 of YLW 15 is con- 
siderably less than the r k > 6 which we estimated 
for this source in Paper II. There may be two 
reasons for this discrepancy. First, we have re- 
examined the older, lower quality data of YLW 15 
and have concluded that CO- band head absorp- 
tion may be present at low confidence in its spec- 
trum (see Figure 2 of Paper II). The presence of 
such absorption would have lowered our estimate 
of the amount of veiling for this source. Second, 
the veiling in the protostar may have decreased 
during the 4 years which elapsed between our two 
observations. Luhman & Rieke (1999) found that 
p Oph sources can undergo large variations in veil- 
ing, and they also derive r k > 1 from their 2 pm 
spectrum of YLW 15. 

We have determined that the rotation velocity 
of the YLW 15 central protostar is v sin i — 50 ± 5 
km s -1 . This was measured by fitting its CO band 
head spectrum (Fig. 3) to spectra of a slowly- 
rotating giant star which had been convolved with 
stellar rotation profiles over a range of velocities. 

The spectra of WL 6 (Figs. 2 and 3) do not 
show any obvious absorption features, so only a 
lower limit to its veiling can be calculated. If the 
central protostar of WL 6 has the same spectral 
type and v sin i rotation velocity as YLW 15, then 
the signal-to-noise of its spectrum (~ 120 in the 
2.2053 - 2.2096 pm region) dictates that WL 6 
must have r k ^ 4.6. This is consistent with the 
lower limit determined in Paper II, r k > 3.8 — 5.5. 
This veiling analysis was based on maximum line 
depths (and not equivalent widths), so we expect 
the veiling of this object to be changed little if it is 
of later spectral type than YLW 15. Our analysis 
of dwarf standards indicates that although equiv- 
alent widths of the 2.2053 - 2.2096 pm features in- 
crease in later types (particularly in the M2 - M6 
range), this is mostly caused by increasing intrin- 
sic line widths (due to pressure broadening), with 
the peak line depths relatively unchanged. How- 
ever, less veiling is required if WL 6 is rotating 
more rapidly than YLW 15, and more is required 
if it is rotating more slowly. Rotation can affect 


5 


the measured line depths and veilings by about a 
factor of 2 over the range v sin i = 0 - 180 km s -1 
(see Paper II). 

4. Discussion 

We have determined that the Class I YSOs WL 
6 and YLW 15 have strong continuum veilings, 
and both have near-IR H 2 emission. Absorption 
lines in the spectrum of YLW 15 have allowed us to 
measure quantitatively its spectral type, r k contin- 
uum veiling, and v sin i rotation velocity. We now 
examine and discuss the physical implications of 
these results given other information known about 
these YSOs. 

4.1. YLW 15 

The SED of YLW 15 rises steeply from the 
near- to the far-IR (2-60 fim) with a bolomet- 
ric luminosity of approximately 10 L 0 (Wilking 
et al. 1989). Its 2-14 fim spectral index (slope 
of its SED) is 1.0 (Bontemps et al. 2001). Andre 
k Montmerle (1994) found that this source also 
has significant A = 1.3 mm spatially-extended flux 
(100 mJy), which they interpreted as arising from 
its circumstellar envelope, confirming its Class I 
nature. Bontemps et al. (1996) found that YLW 
15 has a compact molecular outflow which is seen 
nearly pole-on. The coexistence of molecular out- 
flow and near-IR H 2 emission (Figs. 2 and 3) in 
this object is consistent with the correlation of 
these phenomena found by Greene k Lada (1996). 
GY 263 is a fainter near-IR source (SK = 2.5 mag; 
Barsony et al. 1997) 6" (~ 900 AU projected dis- 
tance) NW of YLW 15. 

However, Girart, Rodriguez, k Curiel (2000) 
recently found that YLW 15 itself is binary in their 
A = 3.6 cm VLA map. 2MASS (second incremen- 
tal release) Point Source Catalog coordinates show 
that the 2 pm source is essentially coincident with 
VLA 2, the weaker, point-like radio source. VLA 
1 is located approximately 076 (~ 100 AU) to the 
NW and is spatially extended (~ 1") SW - NE. 
Girart, Rodriguez, k Curiel (2000) interpret the 
VLA 1 emission to be thermal free-free, arising in 
a protostellar jet. HST / NICMOS near-IR im- 
ages (FI 10W and F160W bands) show that the 
central source is single but embedded in extended 
nebulosity (Allen et al. 2002). Hereafter we as- 
sume that there is a central protostellar engine for 


YLW 15, and it is located at the VLA 2 / near-IR 
position. The VLA 1 source is likely to be emis- 
sion from a jet produced by the central protostar 
and not a (proto-) stellar source itself. 

One of the most intriguing aspects of YLW 15 
is its X-ray emission. This was one of the very few 
Class I protostars which was detected by ROSAT 
(Casanova et al. 1995). Subsequent ROSAT HRI 
(Grosso et al. 1997), ASCA (Tsuboi et al. 2000), 
and Chandra (Imanishi, Koyama, k Tsuboi 2001) 
observations have all shown that its X-ray emis- 
sion has frequent strong flares (~ 10 33 erg s -1 ; see 
Grosso et al. 1997; Imanishi, Koyama, k Tsuboi 
2001). The ASCA and Chandra observations have 
shown that its quiescent and flaring emissions con- 
sist of mostly hard X-rays. The ASCA observation 
captured 3 quasi-periodic flares with a ~ 20 h pe- 
riod which spanned the exposure (Tsuboi et al. 
2000), but the Chandra exposure captured only a 
single flare (~ 3 hr duration) in a 27.8 h exposure. 
Montmerle et al. (2000) have interpreted this flar- 
ing as arising from the interaction of the magnetic 
field of the rapidly rotating protostar (P ~ 20 h, 
near break-up) and its slower rotating circumstel- 
lar disk. 

Our near-IR spectrum of YLW 15 provides fur- 
ther information on the nature of this object. We 
have assigned it a spectral type of K5 IV/ V (see 
§3.2), implying an effective temperature T eff ~ 
4300 K (Tokunaga 2000). Its near-IR spectrum 
(Fig. 3) shows no evidence for binarity; it shows 
no features which could arise in a companion of 
different spectral type. We now derive additional 
stellar and system parameters for this object. 

4- 11. Stellar Parameters 

We used new near-simultaneous near-IR pho- 
tometry (acquired 2001 July 8; see Haisch et al. 
2002) to estimate the extinction towards the cen- 
tral protostar of YLW 15. The observed CIT 
magnitudes of this object were H = 12.53 and 
A = 9.35 (undetected at J) , and we de-reddened 
the H - K color to determine that the extinc- 
tion to the central protostar is A v ~ 40 mag. 
This value is intermediate between the lower limit 
(A v = 34 mag), resulting from de-reddening to 
the maximum CTTS locus color (H - K) 0 = 1.0 
mag (see Meyer, Calvet, k Hillenbrand 1997), and 
the upper limit ( A v = 47 mag), resulting from 
de-reddening to the intrinsic stellar photospheric 
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color (H — K) 0 — 0.11. YLW 15 cannot be ac- 
curately de-reddened to this stellar photospheric 
color because the de-reddened colors of Class I 
YSOs do not intersect stellar loci in near-IR color- 
color diagrams. This is due to the presence of both 
scattered light and infrared excess emission from 
the close-in circumstellar material contained in the 
surrounding accretion disk and infalling envelope 
(see Lada k Adams 1992). Our adopted extinc- 
tion (A v = 40 mag) is consistent with the A v > 31 
mag value derived by Grosso et ah (1997) using 
a similar technique. Our extinction value, along 
with the measured near-IR excess r k = 3.0, the 
IR extinction law, dwarf bolometric corrections, 
and a distance modulus of 6.0 yields an estimated 
stellar luminosity of 3.0 Lq . The 4300 K effec- 
tive temperature of the protostar dictates that its 
stellar radius is 3.1 Rq. Although there are con- 
siderable uncertainties in extinction (and scatter- 
ing), the derived protostellar properties are not 
greatly affected. Varying the A v between its lower 
and upper limits changes the derived stellar lumi- 
nosity by less than a factor of 2 and the derived 
stellar radius by less than 50%. The derived lu- 
minosity and radius are very similar to those pre- 
dicted for a young PMS star of this effective tem- 
perature which is still accreting mass at a high 
rate (M ~ 10” 6 M© yr" 1 ; see Siess, Forestini, k 
Bertout 1999; Tout, Livio, k Bonnell 1999). These 
and other modern PMS models indicate that the 
mass of YLW 15 is approximately 0.5 M©, as de- 
termined from its effective temperature and stel- 
lar luminosity (D’Antona k Mazzitelli 1997; Siess, 
Forestini, k Bertout 1999; Siess, Dufour, k Fores- 
tini 2000). Taken with its derived radius, this dic- 
tates that the surface gravity of YLW 15 is Log 
g ~ 3.1 (cgs units), within the range expected for 
a very young K5 PMS star. 

4-1.2. Accretion 

The derived stellar luminosity is about 30% 
of YLW 15’s bolometric luminosity (Wilking et 
al. 1989), implying that 70% of the total lumi- 
nosity of this protostar is due to mass accretion 
and not radiation of the protostar’s internal en- 
ergy. This accretion luminosity is likely to be 
composed of a component which arises in a cir- 
cumstellar accretion disk as well as a component 
which is generated when accreting matter falls di- 
rectly onto the central protostar. If all matter ac- 


cretes through a circumstellar disk with inner ra- 
dius Rin , then we calculate the accretion luminos- 
ity arising in the disk as GMM/2R in . We assume 
that matter is transported from near the disk edge 
to the stellar surface in magnetospheric accretion 
columns (Koenigl 1991; Shu et al. 1994). We cal- 
culate the accretion luminosity produced by mat- 
ter falling from the disk edge Ri n onto the star as 
GMM/fl* - GMMjRi n . The typical inner hole 
radius of a CTTS circumstellar disk is Ri n ~ 5 R+ 
(Shu et al. 1994; Meyer, Caivet, k Hillenbrand 
1997). Magnetospheric accretion models predict 
that the inner disk radius is closer to the central 
star when matter is accreting at higher rates, such 
as those typical for Class I protostars. DQ Her 
stars are also sometimes interpreted as exhibiting 
this behavior (e.g., see Patterson 1994). Substi- 
tuting the stellar radius (3.1 Rq), stellar mass (0.5 
M©), and mass accretion rate of YLW 15 results 
in a magnetic coupling radius of 3.2 R* for a 1500 
G stellar magnetic field (Shu et al. 1994). 

We derive a mass accretion rate of M = 1.7 x 
10~ 6 M© yr -1 using R in = 3.2 R* and the preced- 
ing terms for the disk and stellar surface accre- 
tion luminosities. This accretion rate is consistent 
with the accretion rate computed for other Class 
I YSOs in this cloud and is within a factor of 2 
of the spherical accretion rate derived from the 
cloud’s gas temperature (M = a 3 /G; see Adams, 
Lada, k Shu 1987). The adopted 0.5 M© stel- 
lar mass of YLW 15 is also consistent with accre- 
tion at this rate over the mean Class I lifetime in 
the cloud (Greene et al. 1994; Luhman k Rieke 
1999). This accretion rate and inner disk radius 
result in accretion luminosities of 1.3 Lq arising 
in the disk and 5.7 Lq arising at the stellar sur- 
face. Thus the stellar surface accretion luminosity 
is nearly twice the luminosity emanating from the 
T eff = 4300 K protostellar photosphere. Caivet 
k Gullbring (1998) have shown that the region 
where an accretion column meets the surface of a 
T Tauri star has an effective temperature of ap- 
proximately 8000 K. From this we calculate that 
the derived 5.7 L© stellar surface accretion lumi- 
nosity arises from an area only 16% of the stellar 
surface. 

4.1.3. Veiling 

Here we investigate the possible sources of ex- 
cess infrared continuum emission and their contri- 
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butions to the overall veiling in the infrared spec- 
trum of YLW 15. We consider three sources for 
the continuum veiling: the accretion luminosity 
generated at the surface of the protostar, the lu- 
minous accretion disk, and the inner infalling en- 
velope. We calculate that the continuum from the 
hot accretion impact sites at the surface of the pro- 
tostar will produce an effective veiling r k = 0.4. 
Although a relatively small fraction of the total 
veiling, this is significantly greater than that pre- 
dicted (r* < 0.05) for accretion impact sites on 
typical class II T Tauri stars whose accretion rates 
are substantially lower (Calvet & Gullbring 1998). 
Accretion models predict that the disk of YLW 15 
{Ldisk = 1.3L 0 ) should contribute an IR veiling 
r k — 1.3 (see Fig. 9 of Greene & Lada 1996). 
The remainder of the object's observed veiling 
(rjt ~ 1.3 of r k =3 total) is likely to arise from 
grains in the circumstellar envelope of the proto- 
star which are warmed by the photons produced 
by the stellar photosphere and at the hot sites of 
mass accretion. This implies that the envelope 
contains a significant amount of dust within ~ 0.3 
AU of the central star (see also Greene & Lada 
1996; Calvet, Hartmann, & Strom 1997). 

We note that the envelope veiling is low enough 
so that any absorption lines originating in the cir- 
cumstellar disk should have been seen in the spec- 
trum of YLW 15. However, the (stellar) dwarf- 
like natures of the YLW 15 absorption features 
imply that they are not formed in a circumstellar 
disk, which would show evidence of much lower 
surface gravity (see Paper I for a discussion and 
examples of FU Ori disk features). We calculate 
that the dust grain destruction distance (where 
T < 1800 K) for YLW 15 is approximately 10 
R* from considering both its protostellar photo- 
spheric and stellar accretion luminosities. This 
implies that the disk is completely gaseous and 
dust-free between R in = 3 R* and 10 /?*. The fact 
that we do not observe absorption lines arising in 
this disk region strongly .implies that the gas in 
this region does not have a sufficient vertical tem- 
perature gradient to form strong spectral lines. 

4-1-4- Rotation 

The v sin i = 50 km s _1 rotation velocity of 
YLW 15 is considerably greater than typical for 
many classical T Tauri stars, about equal to that 
found for flat-spectrum protostars in the p Oph 


cloud (see Paper I), and considerably less than 
the rotation velocity deduced for this object from 
its X-ray emission. If the 20 h X-ray flare period 
observed by Montmerle et al. (2000) was caused 
by rotation as they postulated, then assuming a 
3-1 Rq protostellar radius would imply an equa- 
torial rotation velocity of approximately 180 km 
s" 1 , near break-up. If YLW 15 is rotating this fast, 
then its rotation axis must be only 16° from our 
line of sight. This is consistent with its molecular 
outflow being oriented nearly pole-on (Bontemps 
et al. 1996) if its outflow axis is aligned with its ro- 
tation axis. However, subsequent Chandra X-ray 
observations of YLW 15 do not support the inter- 
pretation that it is rotating with a 20 h period. 
These observations did not confirm flaring with 
a 20 h period; only one flare was detected in a 
10 5 s (27.8 h) observation (Iinanishi, Koyama, & 
Tsuboi 2001). This weakens the possibility that 
the ASCA X-Ray flare period is physically related 
to the protostellar rotation rate. Hereafter we as- 
sume that the rotation axis of YLW 15 is inclined 
57° from our line of sight, the most likely incli- 
nation angle for a random distribution. This im- 
plies that its equatorial rotation velocity is approx- 
imately 60 km s _1 . 

We estimate that the Keplerian co-rotation ra- 
dius of the YLW 15 circumstellar disk is 2.0 R*, 
approximately equal to but slightly less than the 
R = 3 R m magnetic coupling radius. If the cen- 
tral protostar and the disk are sufficiently mag- 
netically coupled in this region, then the disk may 
be regulating the angular momentum (and rota- 
tion velocity) of the star (e.g., Koenigl 1991; Ed- 
wards et al. 1993; Shu et al. 1994). Although the 
observational evidence of such coupling in CTTSs 
has been recently challenged (e.g., see Stassun et 
al. 1999, 2001), the similarity of the co-rotation 
and coupling radii of YLW 15 are consistent with 
stellar spin-down and magnetic star-disk coupling. 
Indeed, we calculate that the magnetic braking 
time scale (see eq. 6 of Hartmann 2002) is only > 
11,000 yr for YLW 15, significantly shorter than 
the - 10 5 yr age of Class I protostars. Thus YLW 
15 is likely to have had more than enough time to 
couple to its disk magnetically. 

Class II YSOs and CTTSs in the p Oph cloud 
(and CTTSs in Taurus; see Edwards et al. 1993) 
are observed to rotate considerably more slowly 
(v sin i ~ 12 km s" 1 ; see Paper I), so if YLW 
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15 is to evolve into a typical Ophiuchi CTTS it 
will likely lose over half of its angular momentum. 
This can result via disk braking if the protostar re- 
mains coupled or recouples to its disk as its mass 
accretion rate drops. If the mass accretion rate of 
YLW 15 suddenly dropped to 10 _7 M© yr _1 , then 
it would re-couple to its disk in only about 50,000 
yr (i.e. eq. 6 of Hartmann 2002). This coupling 
would occur at a radius R ~ 6R, with co-rotation 
speed v sin i — 10 km s -1 (for i = 57°), very sim- 
ilar to the mean rotation velocity observed for p 
Oph Class II YSOs. We note that this 50,000 year 
coupling time is much shorter than the mean age of 
p Oph Class II YSOs (5 - 10 xlO 5 yr), so this can 
explain why Class II YSOs are observed to rotate 
more slowly than Class I and flat-spectrum ones. 
However, more sophisticated modeling is required 
to ascertain the exact nature of this spin-down be- 
havior. 

4.2. WL 6 

Bontemps et al. (2001) found that WL 6 has 
a spectral index of 0.6 (over A = 2 - 14 pm), 
less than that of YLW 15. Wilking et al. (1989) 
also found that the WL 6 spectral index is less 
than YLW 15, but their absolute values and the 
spanned wavelength range are different. WL 6 
has a bolometric luminosity of approximately 2 
L 0 (Walking et al. 1989; Bontemps et al. 2001). 
Neither Andre & Montmerle (1994) or Motte, An- 
dre, & Neri (1998) detected WL 6 at A = 1.3 mm 
(peak F v < 20 mJy in 11" diameter), but both in- 
vestigations regarded it as an embedded protostar 
due to its large spectral index. This source is not 
known to drive a molecular outflow (Bontemps et 
al. 1996), but it does have near-IR H 2 emission 
(Figs. 2 and 3). WL 6 source has no known com- 
panions, and it also appeared single at A = 2 pm 
in the NIRSPEC camera images. 

Kamata et al. (1997) found that WL 6 had 
a hard X-ray spectrum and its flux varied sinu- 
soidally with a 23 h period, which they attributed 
to protostellar rotation. Montmerle et al. (2000) 
also modeled this variation as being due to rota- 
tion, but with a slower period (3.3 d). Imanishi, 
Koyama, & Tsuboi (2001) did not detect this si- 
nusoidal variation, but they did detect two weak 
flares during their 10 5 s Chandra exposure. Taken 
with the results of YLW 15, we interpret these dif- 
ferences to indicate that the X-ray variability of 


WL 6 is not directly related to its protostellar ro- 
tation period. However, this cannot be confirmed 
until the rotation signature of a photospheric fea- 
ture (such as 2.29 pm CO absorption) is observed 
in this object. 

The relatively low luminosity and high veiling 
of WL 6 dictate that its stellar luminosity is very 
low, L , < 0.4Lg, if its effective temperature is 
similar to YLW 15. This is considerably less than 
expected for a young PMS K5 star. Therefore it 
is likely that WL 6 is of later spectral type, lower 
temperature, and lower mass. PMS stars in this 
luminosity range are predicted to have masses ~ 
0.1 Mq (e.g., D’Antona & Mazzitelli 1997; Siess, 
Dufour, & Forestini 2000). Our previous models 
show that it is possible to produce significant veil- 
ing from only a circumstellar accretion disk for a 
low-mass star with this luminosity (see Fig. 9 of 
Greene & Lada 1996). However, it is also plausible 
and possible that matter may also accrete directly 
from the disk to the star, and some veiling could 
also be produced by its circumstellar envelope (all 
as in YLW 15). These limits on the mass, effective 
temperature, and luminosity of WL 6 imply that 
its mass accretion rate is on the order of 10~ 6 A/© 
yr -1 . This rate is consistent with the typical ~ 10° 
yr Class I lifetime and the ~ 0.1 M© mass of WL 
6 . 

5. Summary and Conclusions 

We have presented the first spectrum of a highly 
veiled, strongly accreting Class I protostar (YLW 
15 / IRS 43) which shows late-type photospheric 
absorption features. This spectrum has allowed us 
to diagnose a plethora of details about the stellar, 
accretion, and rotation properties of this proto- 
star and its circumstellar environment. The spec- 
trum of another Class I protostar, WL 6, did not 
show any photospheric absorption features. Both 
of these objects are known to have strong and vari- 
able hard X-ray emission. Our analysis of these 
spectra results in the following conclusions: 

1. We have detected numerous late- type ab- 
sorption features in the photosphere of YLW 15. 
These features include Si (2.0923 pm and 2.2069 
pm), Mg (2.1066 pm), Al (2.1099 pm), Na (2.2062, 
2.2090 pm), Sc (2.2058, 2.2071 pm), Ti (primarily 
2.2217 and 2.2240 pm), and Mg (2.2814 pm) lines 
as well as the CO v = 0 — 2 band. All of these fea- 
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tures are consistent with a single K5 spectral type 
(2"e// = 4300 K) with dwarf-like surface gravity 
(luminosity class IV/ V). 

2. YLW 15 has a low-mass central protostar 
(M 2 ^ 0.5M©) with effective temperature, surface 
gravity, and stellar luminosity similar to that of 
a T Tauri star despite its high bolometric lumi- 
nosity ( Ltoi — 10 Lq). We find that the stellar 
luminosity of YLW 15 is 3 L e and its accretion 
luminosity is 7 £ G . The derived mass accretion 
rate of YLW 15 is M = 1.7 x 10 ~ 6 Mq yr -1 , con- 
sistent with its stellar mass given the lifetimes of 
Class I protostars. 

3. YLW 15 is highly veiled in the near-IR, 
r k > 3. We explore the origin of this veiling and 
evaluate the contributions of the stellar surface ac- 
cretion, circumstellar disk, and circumstellar en- 
velope components. We estimate that the circum- 
stellar disk and envelope contribute equally to the 
veiling each providing ~ 44% of the needed excess 
emission. We find a significant amount (r k ~ 0.4) 
of veiling is provided by the accretion shocks at 
the stellar surface. Although this contribution is 
relatively small (~ 12%) compared to the total, 
it is significantly greater than that produced by 
corresponding processes in typical CTTS. We de- 
termine that emission from the circumstellar en- 
velope is inadequate to veil absorption features 
formed in the inner dust-free circumstellar disk. 
We conclude that there are no strong absorption 
lines formed in the warm inner disk of YLW 15. 

4. The rotation velocity of YLW 15 is v sin 

1 = 50 km s -1 , much faster than that of a typ- 
ical T Tauri star in the p Oph cloud. However, 
this rotation rate does not appear to be correlated 
with the X-ray flarings and variability observed in 
multi-epoch X-ray observations of this source. We 
do find that it is plausible that YLW 15 is magnet- 
ically coupled to its accretion disk, at a radius of 

2 - 3 R . . This suggests - and we verify the plausi- 
bility - that stellar angular momentum decreases 
by about a factor of 2 in only ~ 10^ yr between 
the protostellar and T Tauri evolutionary phases. 

5. The near-IR spectrum of WL 6 does not 
show any obvious absorption features, and we esti- 
mate that the veiling of this YSO is r k >4.6. This 
high veiling and its low luminosity (L(, oi ~ 2 Lq) 
dictate that its central protostar is likely to have 
low stellar luminosity ( L , < 0.4L o ) and to be low 
mass (M ~ 0.1A/ Q ). 
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Fig. 1 — Spectra of MK standards. Three NIR- 
SPEC orders are shown, and strong absorption 
features are identified in each order. Each of the 
three orders contain a number of lines whose ra- 
tios are useful for evaluating the spectral types 
of protostars independent of their veiling. Each 
spectrum has been normalized to a value of 1.0 
and they are offset vertically by a constant in each 
order for separation. 

Fig. 2. — Spectra of Class I protostars YLW 15 
and WL 6. The absorption lines are either weak 
(YLW 15) or undetected (WL 6) due to the large 
continuum veilings (r* > 3) of these objects. Each 
spectrum has been normalized to a value of 1.0 
(over all 3 orders), and W r L 6 has been offset ver- 
tically by 0.2 for in each order for separation. 

Fig. 3. — Spectra of YLW 15 and WL 6 compared 
to HR 5568 (K4 V). The spectra of the protostars 
have been normalized to remove their continuum 
slopes, and 75% of their continuum levels have 
been subtracted to remove the continuum veiling 
amount = 3. Each spectrum has been nor- 
malized to a value of 1.0. The spectra are offset 
vertically by a constant in each order for separa- 
tion. 
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